
ADAPTIVE LINE ENHANCER 

BACKGROUND OF THE INVENTION 
Field Of The Invention 

[0001] The present invention relates to adaptive line enhancers^ 
and to methods for adaptive line enhancement. Applications for the 
5 invention lie in the fields of radar, sonar, communications and 
other related disciplines where digital signal processing may be 
required. 

BACKCROUND OF THE INVENTIO N Description Of The Related Art 
10 [0002] Detection of sinusoidal signals immersed in noise is a 
fundamental problem in signal processing. The retrieval of 
sinusoidal or other narrow-band signals which may have been 
significantly attenuated, frequency shifted because of Doppler 
effects and corrupted by interference and noise^ has conventionally 
15 been carried out using analysis of the signal in the frequency 

domain. This requires the input signal to be Fourier transformed. 
Once the signal has been Fourier transformed, the strongest 
spectral component can be detected and a filter designed to either 
enhance or reject this frequency. For the detection of sinusoids 
20 with time-varying frequencies, a Fourier transform with sliding 
windows can be used. Despite the availability of algorithms^ such 
as the fast Fourier transform (FFT) which are computationally 
efficient when compared to a direct implementation of the discrete 
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Fourier transform, the frequency domain analysis of the input 
signal is still relatively inefficient when compared to adaptive 
line enhancement techniques. 

[0003] Adaptive line enhancement is an alternative technique to 
5 frequency domain analysis based on FFTs. It has been shown (B. 
Widrow and S.D. Stearns, u Adaptive Signal Processing" , Prentice- 
Hall, 1985) that Adaptive Line Enhancers (ALEs) require fewer 
computations than FFT FFT- based techniques and^_ in certain 
circumstances^, can be more sensitive detectors of sinusoids. The 

10 ALE consists of a filter, and an adaptation rule for changing some 
feature of the filter's frequency response characteristics. Various 
combinations of filters and adaptation rules have been proposed, 
with the most recently reported embodiments comprising a lattice 
Gray-Markel adaptive notch filter and adaptation rules based on a 

15 simplified gradient technique (N.I. Cho, C.-H. Choi and S.U. Lee, 
"Adaptive Line Enhancement by Using an IIR Lattice Notch Filter," 
IEEE Trans. Acoust . , Speech, Signal Processing, vol. 37, Apr. 1989; 
P. A. Regalia, w An Improved Lattice-Based Adaptive IIR Notch 
filter," IEEE Trans. Signal Processing, vol. 39, pp. 2124-2128, 

20 Sept. 1991). It has been shown that such ALEs provide better 

convergence to the frequency of interest than previous designs and^ 
in addition^ are less sensitive to the finite word length effects 
which occur in any digital processor. 

[0004] The transfer function of the Gray-Markel lattice notch 
25 filter is expressed as: 
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rr N(z) A + a. \ + 2k^+z- 2 ^ x 
TJ, . =— — = ( ) - : r- (Equation A) 

where k 0 determines the notch frequency and where a determines the 
5 bandwidth. The notch frequency determining variable k 0 should 
converge to -cos(6; 0 ) to reject a sinusoid with frequency w 0 . This 
filter has zeros on the unit circle at z 0 =e ±0) ° , where to 0 = cos _1 (-A: 0 ) . 
The -3dB attenuation bandwidth bandwidth BW of the magnitude 
response of the Gray-Markel lattice notch filter is determined by 
10 the following equation: 

BW = cos~ ] (-^- T ) 
\ + a 

1 + a 

[0005] A slight gain correction of ( ) is needed to achieve 

15 unity gain in the passband. 

[0006] The bandwidth and the notch frequency can be controlled 
separately by changing k 0 and a. This filter structure can easily 
be implemented using either a direct form realization or a lattice 
filter structure based on wave digital filters (WDFs) (A. Fettweis 
20 and H. Levin and A. Sedlmeyer, "Wave Digital Lattice Filters," Int. 
J. Circuit Theory Applicat, vol. 2, no. 2, pp. 203-211, June 1974; 
A. Fettweis, "Wave Digital Lattice Filters: Theory and practice 

(invited paper) , " Proc. IEEE, vol. 74, pp. 270-327, Feb. 1986). 
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[0007] Referring to Fig. 1, there is shown a block diagram 
showing functional elements for implementing a known wave digital 
filter realization of the Gray-Markel notch filter response, the 
transfer function of which is given in Equation A. 
5 [0008] In Fig. lj_ there is shown an input 110, a first dynamic 
adapter block 120, a second dynamic adapter block 130, a summing 
block 140, an amplifier block 150, an output 160, a notch bandwidth 
determining block 170 and a notch frequency determining block 180. 

10 [0009] An input signal 110 is fed to a first input of the 

summing block 140 and to a first input of the first dynamic adapter 
block 120. A first output of the first dynamic adapter block 120 is 
fed to a second input of the summing block 140. The output of the 
summing block 140, comprising the result of the addition of input 

15 110 and a first output of the first dynamic adapter block 120^ is 
fed to the input of the amplifier block 150. Amplifier The 
amplifier block 150 has a fixed amplitude gain of 0.5. This gain is 
achieved by a bit -shift operation^ and thuSj_ does not require a 
multiplier. The output of the amplifier block 150 becomes the 

20 output signal 160. 

[0010] A second output of the first dynamic adapter block 120 is 
left unconnected. A third output of the first dynamic adapter block 
120 is fed to a first input of the second dynamic adapter block 
13 0, said third outputs in combination wit h the second dynamic 
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adapter adapter 130^ forming a feedback path aroun d the first 
dynamic adapter block 120. 

[0011] A first output of the second dynamic adapter block 130 is 
fed back to a second input of the first dynamic adapter block 120. 
5 The output of the notch bandwidth determining block 170 is fed to a 
third input of the first dynamic adapter block 120. 
[0012] A second output of the second dynamic adapter block 13 0 
is left unconnected. A third output of the second dynamic adapter 
block 130 is fed back to a second input of the second dynamic 
10 adapter block 130. The output of the notch frequency determining 
block 180 is fed to a third input of the second dynamic adapter 
block 13 0. 

[0013] Referring to Fig. 2 , there io ohown shows a block diagram 
15 showing functional elements for implementing the dynamic adapters 
120 and 130. 

[0014] In Fig. 2j_ there is shown a first input 210, a second 
input 220, a third input 230, a first subtracter block 240, a 
multiplier 250, a second subtracter block 260, a third subtracter 
20 block 270, a delay block 280, a first output 285, a second output 
290 and a third output 295. 

[0015] The first input 210 is fed to the positive input terminal 
of the first subtracter block 240 and to the negative input 
terminal of the third subtracter block 270. A second input 220 is 
2 5 fed to the negative input terminal e #of the first subtracter block 
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240 and the negative input terminal of the second subtracter block 
260. The output of the first subtracter block 240, comprising the 
difference of the first input 210 an d the second input 220^_ is fed 
to a first input terminal of the multiplier 250, A third input 230 
5 is fed to a second input terminal of the multiplier 250. The output 
of the multiplier 250, comprising the product of the third input 
230 and the output of the first subtracter block 240, is fed to the 
positive input terminal of the second subtracter block 260 and to 
the positive input terminal of the third subtracter block 270. The 

10 output of the second subtracter block 260, comprising the 

difference of the output of the multiplier block 250 and the second 
input 220, becomes the first output 285. The output of the third 
subtracter block 270, comprising the difference of the output of 
the multiplier 260 and the first input 210j_ becomes the second 

15 output 290 and is also fed to the delay block 280. Delay The delay 
block 280 delays the signal by one sample period then feeds it to 
the third output 295. 

[0016] The time domain response of the dynamic adapter can be 
evaluated. From an initial state where the values of all inputs and 

20 outputs are zero, a train of pulses a 1# a 2 ,...,a n , a n+1 ,... is applied to 
the first input 210, a train of pulses b 1# b 2 ,...,b n/ b n+1 ,... is applied 
to the second input 220^ and the constant K is applied to the third 
input 230 then . Then the outputs at the n th time step become: 
First output 285 (a^bj X K - b n 

25 Second output 290 (a^b^ X K - ^ 
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Third output 295 (a^-b^) X K - 

[0017] The time domain response of the dynamic adapter^ shown in 
feedback configuration in Fig. 5 (and in the dashed box in Fig. 1)^ 
5 can be evaluated. From an initial state where the values of all 
inputs and outputs are zero, when a train of pulses u(l), 
u(2) ,...,u(n) , u(n+l) # ... is applied to the input 510 and the constant 
K (520) is applied to third input of dynamic adapter block 530 , the 
equation for the output 540 at the n th time step becomes: 

10 

Output 540 y(n) = Kxu(n) + u(n-\)-Kxy(n-\) 

[0018] This corresponds to an all -pass transfer function in the 
z- transform domain of 

15 

rr = *±£? 

±1 adaptor j _j_ 

[0019] It is desirable to reduce the computational complexity of 
ALEs. 

20 

SUMMARY OF THE INVENTION 
[0020] It is an object of the present invention to provide an 
apparatus for adaptive line enhancement using a computationally 
efficient adaptation mechanism for an adaptive Gray-Markel lattice 
25 notch filter. 
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[0021] It is a further object of the present invention to 
provide a method for adaptive line enhancement using a 
computationally efficient adaptation mechanism for an adaptive 
5 Gray-Markel lattice notch filter. 

[0022] The inventor has determined that an ALE of reduced 
computational complexity can be achieved by applying a new 
adaptation algorithm. An algorithm for adapting the notch frequency 
determining variable, k 0/ by the recursive Gauss -Newton algorithm 
10 with simplified gradient and sgn-sgn adaptation rule is as follows: 
Algorithm 1: 
Input signal is U(z) 

Gray-Markel notch filter: ff, . = — — = ( ) ; T 



15 



Output signal is ^( z )~Hi a ai Ce 



Attempt to minimize the expected value of the energy of the output 
E(Y 2 ) 

Filter Operations (for each input sample) : 

Calculate notch filter output y(n) using WDF filter structure (see 
Fig. 1) 
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Adapt k according to Equation C (below) in which x(n-l) is the 
simplified gradient of the current sample and relates to input u(n) 
as shown in Equation B. 

5 x(n) = — — u(n) = u(n) - k(n)(\ + a)x(n - 1) - ax{n - 2) 
D(z) 

k(n + 1) = k(n) - sgn[x(n - l)y(ri)]/u 
p. is the adaptation constant . 

10 

Stability monitoring: clip k(n+l) within range ] -1 1[ 

[0023] The adaptation constant p determines the rate of 
convergence of the algorithm on k 0/ and also puts bounds upon the 
15 achievable mean accuracy of the estimation of kO . 

[0024] To avoid the calculation of the simplified gradient x(n- 
1) needed for the update of k 0 , the internal variables of the WDF 
filter structure of the Gray-Markel notch filter are studied. 
[0025] The table below shows a z- transform domain transfer 
20 function of the internal variables in the signal flow graph of the 
wave digital filter structure shown in Fig. 1. 



(Equation B) 



(Equation C) 



Internal Variable Transfer function 

Input 1 

25 
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Output 

Outl Adapter— 1 
5 Out 3 Adapt er-1 = Inl Adapt er-2 

In2 Adapt er-1 = Outl Adapt er-2 

In2 Adaptor 2 = Out 3 Adapt er-2 

< 

10 Out 2 Adapter-2 



2 'l + k 0 (l + a)z x +az 2 

1+g a + k 0 (l+a)z~ l + z~ 2 
~2~ 1 + k 0 (1 + a)z' x + az~ 2 



(a-l)z 



(a-l)z 



-i 



l + fc 0 z~ 



\ + k 0 (\+a)z 1 +az 2 



k 0 +z 1 



(g-l)fe^fl^ " t+flg " 2 
l + A: 0 (l + a)z _1 +az" 2 

lH-A: 0 (l + a)z" 1 +oz" 2 



[0026] 



It can be seen from the definition of simplified gradient 



x(n-l) given in Equation B and from the z-transform domain transfer 
functions of the internal variables of the wave digital filter 
15 shown in Fig. 1^ that the third output of second dynamic adapter 
block 13 0 corresponds to (a-1) (k 0 -l)x(n-2) j_ and that the second 
output of second dynamic adapter block 130 corresponds to (a-1) (k 0 - 
l)x(n-l) . 

[0027] Since both |a|<l and |k|<l as a fundamental requirement 
20 for stability of the Gray-Markel filter, the product (a-1) (k 0 -l) is 
always positive. Therefore: 



sgn [x(n-l)y(n) ] =sgn [y (n) ] sgn [Out 2 Adaptor 2 ] 



(Equation D) 
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_This equation makes the calculation of simplified 



gradients unnecessary and leads to a new algorithm for an ALE using 
a low complexity adaptive lattice notch filter based on wave 
digital filters. The algorithm is as follows: 

5 

Algorithm for Adaptive Line Enhancer using a Low Complexity 
Adaptive Lattice Notch Filter based on Wave Digital Filters: 

Input signal is u(n) with n starting at time 0 (U(z) in frequency 
10 domain notation) 

N(z) ,l + a x l + 2£ 0 z _1 + z~ 2 
Gray-Markel notch filter: WV =— — = ( ) : T 

Output signal is Y(z) = fj^ x U(z) 

15 

Attempt to minimize the expected value of the energy of the output 
E(Y 2 ) 

Initialization: 
20 Initialize Out^Adaptor x {-\) 
Ini t ial i ze Out l Adaptor 2 (-1) 

Filter Operations (for each input sample u(n)): 
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Calculate notch filter output y(n) using WDF filter structure with 
input u(n) (see Fig. 1) : 

Out 3 Adaptor x (n) = Out 2 Adaptor x (n - 1) 
5 Out 3 Adaptor 2 (n) = Out 2 Adaptor 2 {n - 1) 

Out x Adaptor 2 (n) = k(n) x [Out 3 Adaptor x (n) - Out 3 Adaptor 2 («)]- Out 3 Adaptor 2 (n) 
Out 2 Adaptor 2 (n) = k(n) x [Out 3 Adaptor x (n) - Out 3 Adaptor 2 (n)]- Out 3 Adaptor x (n) 
Out x Adaptor x (n) = ax [u(n) - Out x Adaptor 2 («)]- Out x Adaptor 2 (n) 
Out 2 Adaptor x (n) = ax [u(n) - Out x Adaptor 2 (n)]-u(n) 
1 0 y(n) = 0.5 x [u(n) + Out x Adaptor x («)]- u(n) 

Update of the variable k determining the notch frequency (see Fig. 
3) : 

k(n + 1) = k{n) - sgn[y(n)]sgp[Out 2 Adaptor 2 ]x ju 

15 

p is the adaptation constant. 

Thus Out 2 Adaptor 2 can be used as an update function (UPDATEFN) for 
an update of the variable k used to determine the notch frequency 
20 in a given sampling period. As referred to above, UPDATEFN 

therefore has a transfer function in the z-transform domain, in 
this case x in an n th sampling period of: 

(a -!)(£(*) -IK 1 
1 + £(«)(! + a)z _1 + cxz 2 
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Stability monitoring: clip k(n+l) within range ]-l 1 [ . 

[0028] According to the present invention^ in a first aspect, 
there is provided an adaptive line enhancer comprising an adaptive 
5 Gray-Markel lattice notch filter having an adaptive notch 

frequency, in which the notch frequency is determined according to 
a notch frequency variable k, characterized in that the value of k 
for the n+l th sample period is determined according to the following 
equation: 

10 

k(n + 1) = k(n) - sgnbOOl s&i[UPDATEFN]x M 

in which y(n) is the notch filter output, p is the adaptation 
constant, and UPDATEFN has a transfer function in the z-transform 
15 domain of: 

(a-l)(k(n)-l)z- 1 
l + k(ri)(l + a)z~ l +az~ 2 

in which a determines the bandwidth and k(n) is a variable for 

20 determining the current notch frequency. 

[0029] According to the present invention^ in a second aspect, 
there is provided a method for adaptive line enhancement, 
comprising adaptive line enhancing an adaptive Gray-Markel lattice 
notch filter with an adaptive notch frequency, in which the notch 

25 frequency is determined according to a notch frequency variable k, 
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characterized in that the value of k for the n+1 sample period is 
determined according to the following equation: 

k(n + 1) = k(n) - s&i[y(n)]s&i[UPDATEFN]x M 

5 

in which y(n) is the notch filter output, p. is the adaptation 
constant, and UPDATEFN has a transfer function in the z- transform 
domain of: 

iQ (a-l)(k(n)-l)z- 1 
1 + k(n)(l + a)z~ l +az~ 2 

in which a determines the bandwidth and k(n) determines the current 
notch frequency. 

15 [0030] This algorithm for adapting the notch frequency enables 

UPDATEFN and the notch frequency variable to be directly calculated 
from knowledge of internal variables of the wave digital filter 
[0031] The update rule for k(n) in this algorithm is very simple 
and requires, when compared with the prior art notch filter of Fig. 

20 1, only an extra computational load of one addition, two sign 

operators and an EXOR operator (for doing the multiplication of the 
two signs) . Together with the operation for calculating the WDF 
filter structure^ this makes a total of two multiplications, eight 
additions, one bit bit- shift, two sign operators and one EXOR 

25 operator^ for each processed input sample. 
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[0032] This represents a saving of three multiplications and two 
additions which would otherwise be needed to calculate the 
simplified gradient each sample period, according to Equation B. 

Claims 2 5 define advantageous apparatus for putting the 

5 present invention in to effect, 

Claims 7 10 define advantageous ways in which the method 

of the present invention may be implemented. 

BRIEF DESCRIPTION OF THE DRAWINGS 
10 [0033] For a better understanding of the invention, and to show 
how an embodiment of the same may be carried into effect, reference 
will now be made, by way of example, to the accompanying 
diagrammatic drawings in which: 

[0034] Fig. 1 is a schematic drawing b lock diagram illustrating 
15 means of implementing a prior art Gray-Markel notch filter using a 
wave digital filter; 

[0035] Fig. 2 shows a schematic drawing block diagram 
illustrating means of implementing the "dynamic adapter" shown in 
Fig. 1; 

20 [0036] Fig. 3 shows a schematic drawing block diagram 

illustrating an adaptive line enhancer according to the present 
invention; 

[0037] Fig. 4 shows the results of a frequency hop experiment 
used to evaluate the effectiveness of an embodiment of the present 
25 invention; and 
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[0038] Fig. 5 shows a schematic drawing b lock diagram 
illustrating the "dynamic adapter" in a feedback configuration. 

DESCRIPTION OF THE PREFERRED ENBOD IMENTS 
5 [0039] Referring now to Fig. 3, there is shown a block diagram 
that is identical to the Fig. 1 schematic block diagram except that 
jrfc — the block diagram of Fig. 3 includes additional feedback 
elements that implement for implementing the low complexity 
adaptation algorithm that forms the basis of the invention. 

10 [0040] In Fig. 3^_ there is shown an input 305 u(n), a first 
dynamic adapter block 310, a second dynamic adapter block 315, a 
first summing block 320, an amplifier block 325, an output 330, a 
notch bandwidth determining block 335, a notch frequency 
determining block 340, a first signum function block 345, a second 

15 signum function block 350, a first multiplier 355, a second 

multiplier 360, an adaptation speed determining block 365, a second 
summing block 370, an amplitude limiting block 375 and a delay 
block 380. 

[0041] An input signal 305 is fed through a first input of the 
20 first summing block 320 and to a first input of the first dynamic 
adapter block 310. A first output of the first dynamic adapter 
block 310 is fed to a second input of the first summing block 320. 
The output of the first summing block 320, comprising the result of 
the addition of input 3 05 and a first output of the first dynamic 
25 adapter block 310j_ is fed to be— an input of the amplifier block 
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325. Amplifier The amplifier block 325 has a fixed amplitude gain 
of 0.5. This gain is achieved by a bit-shift operation^, and thus^ 
does not require a multiplier. The output of the amplifier block 
325 becomes the output signal 330^ and is also fed to the input of 
5 second signum function block 350. A second output of the first 
dynamic adapter block 310 is left unconnected. A third output of 
the first dynamic adapter block 310 is fed to a first input of the 
second dynamic adapter block 315. A first output of the second 
dynamic adapter block 315 is fed back to a second input of the 

10 first dynamic adapter block 310. The output of the notch bandwidth 
determining block 335 is fed to a third input of the first dynamic 
adapter block 310. A second output (Out2) of the second dynamic 
adapter block 315 is fed to the input of the first signum function 
block 345. A third output of the second dynamic adapter block 315 

15 is fed back to a second input of the second dynamic adapter block 
315. The output of the first signum function block 345 is fed to a 
first input of the first multiplier 355, and the output of the 
second signum function block 350 is fed to a second input of the 
first multiplier 355. The output of the first multiplier 355, 

20 comprising the product of the output of the first signum function 
block 345 and the second signum function block 350^ is fed to a 
first input of the second multiplier 360. The output of an-the 
adaptation speed determining block 365 is fed to a second input of 
the second multiplier 360. The output of the second multiplier 360, 

25 comprising the product of the output of the first multiplier block 



N: \UserPublic\GO\SX10GOA0 .RED 17 



* 



PHBE 000020 

355 and the output of the adaptation speed determining block 365j_ 
is fed to a first input of the second summing block 370, The output 
of the notch frequency determining block 340 is fed to a third 
input of the second dynamic adapter block 315 and to a second input 
5 of the second summing block 370. The output of the second summing 
block 370, comprising the sum of the output of the second 
multiplier 360 and the output of the notch frequency determining 
block 340^_ is fed to the input of the amplitude (saturation) 
limiting block 375. The output of the amplitude limiting block 375 
10 is fed to the input of the delay block 380. The output of the delay 
block 380 becomes the updated value of the notch frequency 
determining block 34 0 and^_ accordingly^ is fed to a third input of 
the second dynamic adapter block 315. 

[0042] Amplitude limiting block 375 prevents k(n+l) from 
15 becoming M or ^ -1. When |k(n+l) the notch filter becomes 

unstable. To prevent instability^ k(n+l) is clipped in— to the open 

interval ] -1 1[. This is done as follows: 

If k(n+l) ^ClipValue then k (n+1) =clipvalue 
If k(n+l)^ ClipValue then k (n+1 ) =-clipvalue 

20 

[0043] Wit h where clipvalue being is slightly less than 1, e.g.j_ 

0.999. This is also referred to as stability monitoring. 
[0044] The second output Out 2 of the second dynamic adapter 315 
is used to generate the k(n+l) value^ used as a variable to 
25 determine the update for the adaptive coefficient determining the 
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notch frequency. The signum of 0utput2 of the second dynamic 
adapter 315 is generated b y the first signum block 345, which is 
multiplied b y the first multiplier 355 with the signum of the 
output of the amplifier block 325 (which is y(n) , the notch filter 
5 output) , the signum of the output of the amplifier block 325 being 
carried out b y the second signum block 350. This^_ therefore^ 
generates sgn [y (n) ] sgn [Out2] as the output of first multiplier 355, 
This is then multiplied by adaptation constant p. fro m the 
adaptation spee d determining block 365 at the second multiplier 360 
10 and subtracted from the current k(n) to generate k(n+l) . 

[0045] ThuSj_ the second output Out2 of the second dynamic 
adapter 315 is used as an update function (UPDATEFN) . As shown in 
the table above (reference Output 2 Adaptor 2 ) j_ UPDATEFN has a transfer 
15 function in the z-transfer domain, for an n th sample of: 

(a-l)(k(n)-l)z- 1 
l + k(n)(l + a)z- x +az~ 2 

[0046] The embodiment of the present invention^ represented 
20 diagrammatically by the block diagram shown in Fig. 3j_ has 
significant advantages over prior realizations of ALEs, 
particularly in terms of minimizing the amount of hardware needed 
to carry out the ALE procedure and minimizing the computational 
load needed to carry out the ALE procedure on any digital 
25 processor. Fig. 4 shows the results of a frequency hop experiment 
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in which the input signal^ supplied to the embodiment of the 
invention shown in Fig. 3j_ is a sine wave immersed in white noise 
and sampled at a sampling rate f s equal to 16kHz. The frequency of 
the sine wave changes randomly every 1000 samples. The first graph 
5 (from top to bottom) shows the desired frequency (des freq) with 
oc=0.8, p=0.005 and SNR=23dB. The second graph shows the estimated 
frequency (est freq) using an ALE of the described embodiment of 
the present invention. The third graph corresponds to the first 
graph except that a=0.7, p=0.001 and SNR=4.9dB. The fourth graph 
10 shows the corresponding estimated frequency using an ALE of the 
described embodiment of the present invention. 

The embodiment of the present invention is used to make 

an estimate of this desired frequency /^. esftm .( w ) for each time step 
by using: 

15 

ffr^estonW = ^cos" 1 ^*)] (Equation E) 

[0047] To quantify the achievable accuracy of the frequency 
estimation given by this algorithm in case of the frequency hop 
20 experiment, it can be proved p roven that the standard deviation of 
the estimated frequency a> 0 = cos -1 (-£ 0 ) shows following proportional 
relationship: 



Vfret.estim. « / , (Equation F) 



real 

25 
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[0048] The estimate of the desired frequency will be converged 
on— to the desired frequency in a time showing the following 
proportional relationship: 

5 T conv « — samples (Equation G) 

M 

[0049] By choosing appropriate values for a and p by inspection 
of equations F and G and/or basic experimentation, useful adaptive 
sinusoid tracking can be achieved for various oignal signal -fe e — to- 

10 noise ratios^ as is demonstrated in Fig. 4. 

Accordingly^ the variable k^ used for determining the 
notch frequency^ is itself determined by an output of the second 
dynamic adapter i.e.j_ internally rather than externally. 
[0050] All of the features disclosed in this specification 

15 (including any accompanying claims, abstract and drawings) , and/or 
all of the steps of any method or process so disclosed, may be 
combined in any combination, except combinations where at least 
some of such features and/or steps are mutually exclusive, 
[0051] Each feature disclosed in this specification (including 

20 any accompanying claims, abstract and drawings) , may be replaced by 
alternative features serving the same, equivalent or similar 
purpose, unless expressly stated otherwise. Thus, unless expressly 
stated otherwise, each feature disclosed is one example only of a 
generic series of equivalent or similar features. 
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[0052] The invention is not restricted to the details of the 
foregoing embodiment ( s ) . The invention extend to any novel one, ■ 
any novel combination, of the features disclosed in this 
specification (including any accompanying claims, abstract and 
drawings) , or to any novel one, or any novel combination, of the 
steps of any method or process so disclosed. 
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ABSTRACT- OF THE DISCLOSURE 

The prcocnt invention dioclooco aAn improved adaptive 
line enhancer comprising includes an adaptive Gray-Markel lattice 
5 notch filter having an adaptive notch frequency, in which the notch 
frequency is determined according to a notch frequency variable k- 
charac tori zed in that t . T he value of k for the n+l th sample period 
is determined according to the following equation: 
k(n + 1) = k(n) - s&i[y(n)]sgn[UPDATEFN]x M 
10 in which y(n) is the notch filter output, p is the adaptation 

constant, and UPDATE FN has a transfer function in the z- transform 
domain of : 

l + k(n)(l + a)z~ l +az~ 2 

15 

in which a determines the bandwidth and k(n) is a variable for 
determining the current notch frequency, 

A corresponding method io also diocloocd. 

__The algorithm for adapting the notch frequency enables 

20 the notch frequency to be directly calculated from knowledge of 
internal variables of the wave digital filter. 
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